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Abstract 

This paper concerns the 3-dimensional Lagrangian Navier-Stokes a model and the lim- 
iting Navier-Stokes system on smooth bounded domains with a class of vorticity-slip 
boundary conditions and the Navier-slip boundary conditions. It establishes the spec- 
trum properties and regularity estimates of the associated Stokes operators, the local 
well-posedness of the strong solution and global existence of weak solutions for initial 
boundary value problems for such systems. Furthermore, the vanishing a limit to a 
weak solution of the corresponding initial-boundary value problem of the Navier-Stokes 
system is proved and a rate of convergence is shown for the strong solution. 
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1 Introduction 



The Lagragian Navier-Stokes a model (LNS-a) as a regularization system of the Navier- 
Stokes equations (NS) is given by 

d t v - Av + T a v ■ Vv + V(T a v) T -v + Vp = (1.1) 
V-v = (1.2) 

which describes large scale fluid motions in the turbulence theory, where T a v = u is a 
filtered version of the velocity v determined usually by 

u — aAu = v (1-3) 
V-u = (1.4) 

with a > being a constant. This filter u is also called the averaged velocity. The 
system can be regarded as a system for this filter, and is also called the Lagrangian 
averaged Navier-Stokes equations (LANS). The ideal case, called the Lagrangian aver- 
aged Euler equations (LAE) or Camass-Holm equations, was first introduced in [144127] . 
The viscosity was added in [15\ \W[ [28] yielding the LANS which is sometime called 
viscous Camass-Holm equations. 

The global well-posedness for the LANS was first obtained in [21] for periodic bound- 
ary conditions. The convergence of its solutions to that of the NS equations and the 
continuity of attractors when a — > are also considered there. 

For bounded domains, the situation becomes more complicated since the LANS is a 
4th odder system for the filter u, and only the no-slip boundary condition u = on the 
boundary was considered by [37] under the assumption that Au = —PAu = on the 
boundary with P being the Leray projection operator. The boundary effects related 
to such a boundary condition were analyzed in [29] . We also refer [211 12H ESI E7J for 
more details along this line. 

On the other hand, the LNS-a model emphasizes the system (1.1)-(1.4) as equations 
for the physical velocity v, which is a regularized system of the NS equations by filtering 
some part of the nonlinearity through a global quantity which is then called filtered 
velocity (see [24] and the references therein). There are many filtered formulations, 
which thus lead to many a models(see [12} 125] for instance). It is also mentioned in [18] 
in the stochastic Lagrangian derivation of (1.1), (1-2) that any translation-invariant 
filter u = T a v may be adaptable. 

Although, there is no any serious difference between the two aspects for the equations 
(1.1), (1.2) filtered by (1.3), (1.4) in domains without boundary, the situation may be 
different for domains with boundaries. To our knowledge, very little is known to the 
LNS-a models in domains with boundaries from this point of view. 

In this paper, we investigate the initial boundary value problem for the LNS-a model 
(1.1), (1.2) in the following equivalent form 

d t v - Av + V x v x T a v + Vp = in (1.5) 
V-v = 0mn (1.6) 







in a smooth bounded domain with the property that both f2 and <9f2 have only finite 
many simply connected components, where V- and Vx denote the div and curl oper- 
ator, respectively. 

Once the filter mapping T a is given, equations (1.5) and (1.6) become a Navier- 
Stokes type system for v, and for which, some boundary conditions are needed. Here 
we consider the following vorticity-slip boundary condition (VSB): 

v ■ n = 0, n x V x v = (3v on <9fi (1.7) 

Since there is a boundary, the filter u = T a v can not be determined by solving (1.3) and 
(1.4). Some boundary conditions are also needed. We propose that the filter u = T a v 
be determined by solving the following Stokes boundary value problem 

u — aAu + Vp = v in Q (1-8) 
V-u = 0inO (1.9) 

with the VSB: 

u ■ n = 0, n x V x u = j3u on dfl (1-10) 

We also consider the associated boundary value problem for the Navier-Stokes equations 

d t v° - uAv° + (V x v°) x v° + Vp = (1.11) 
V-v° = (1.12) 

with the corresponding boundary conditions (1.7) as a singular limit problem by pass- 
ing to the vanishing a limit in (1.5)-(1.10). 

The equivalence between (1.1) and (1.5) follows from the vector formula 

V(u-v) = u-Vv + Vu T -v-V xvxu (1.13) 
for any divergence free vectors u and v. 

There have been extensive studies of the Navier-Stokes systems on bounded domains 
with various boundary conditions, such as the well known no-slip condition and various 
slip boundary conditions. In particular, substantial understanding has been achieved 
for the well-posedness of initial boundary value problems for the Navier-Stokes system 
with these boundary conditions and problems of vanishing viscosity limit and boundary 
layers, see [3 H CDS EH [201 E3 [331 G21 [351 H3 HE] and the references therein. Note 
that the no-slip boundary condition corresponds to our VSB with f3 = oo. Yet one of 
the main motivations for the proposed VSB is its relation to the well known Navier-slip 
boundary condition (see [H El [H EJ [30], [331 ED HH] and the references therein). Indeed, 
the Navier-slip boundary condition (NSB) says that the fluid at the boundary is allowed 
to slip and the slip velocity is proportional to the shear stress (see [39]), i.e., 

v ■ n = 0, 2((S(v)n)) T = - 7 t; T on dtl (1.14) 

where 2S(v) = (Vi> + (Vi> ) T ) is the stress tensor. Note that 

(2(S(v)n) -(Vx»)x n) T = GD{v) T on dtt (1.15) 
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where GD{v) = —2S(n)v is the lower order term due to the geometry of the bound- 
ary, see lemma 3.10. below. In the special case that the boundary <90 is flat, one has 
GD(v) = 0. Thus the VSB (1.7) coincides with NSB (1.14). It should be mentioned 
that as far as we know, all the previous physical and numerical studies concerning the 
NSB deal with only the case of flat boundaries [TJ [HJ I3T1 [32l l4"2"j [45] . Another main 
motivation for the proposed VSB (1.7) and (1.10) is that the vorticity formulations of 
the fluid equations have played important roles in analyzing fluid motions, and suitable 
boundary conditions on the vorticity should be important for such formulations, see 
[H 13 CG2 HH [36] and the references therein. For example, the equivalent vorticity form 
of the NSB conditions are crucial in the studies of the corresponding boundary value 
problems in |10[ 117] , and the VSB (with /3 = 0) was found very useful to understand 
the vanishing viscosity limit problem of the Navier-Stokes equations in [31 [35| |46| 147] . 
It is hoped that the VSB conditions proposed here can share light on understanding 
the fluid motions in bounded domains. 

The rest of the paper is organized as follows: First, as a preparation, we present 
in the next section a 1? version of the general Hodge decomposition theory that was 
stated in [TlJ for smooth vector spaces, which will be used to study the Stokes problems 
associated with various slip boundary conditions. Then we give general and systematic 
results on well-posedness and spectrum properties of the Stokes operators associated 
with various VSB and NSB conditions in section 3. Our results apply to domains 
with general topology. It should be mentioned that all the previous analysis deals 
with only the NSB conditions in some special cases. Based on the properties of the 
Stokes operators, in section 4, we can formulate the initial boundary value problem of 
the LNS-a model, (1.5)-(1.10), together with the limit problem of the NS equations, 
(1.11), (1.12), (1.7), to be a series of abstract equations in a Hilbert space for the pa- 
rameter a G [0, oo). In section 5, we study the well-posedness of the weak solutions for 
the LNS-a equations with the VSB conditions for each a > 0, by the Galerkin method. 
The local well-posedness, theorem 5.1., is obtained by direct estimates on the velocity 
v, while the global theory, theorem 5.2., is proved by combining energy estimates on 
both the velocity field v and the filter u. Note that our approach is somewhat dif- 
ferent from [21\ [37] in emphasizing the velocity v but not the filter u. In section 6, 
we investigate the vanishing a limit of solutions of the initial boundary value problem 
of the LNS-a equations with VSB condition to the corresponding solutions of the NS 
equations. The global in time convergence of weak solutions is obtained in theorem 
6.1. similar to periodic case in |21| . while local in time convergence of strong solutions 
is given in theorem 6.2.. The existence of the global weak solutions and local unique 
strong solution for the NS equations with corresponding VSB condition are then fol- 
lowed. Furthermore, some estimates on convergence rates are given in theorem 6.3.. 
Finally, we present some generalizations in section 7. In particular, a parallel theory 
holds for the NSB condition. 



2 Preparations and Hodge decompositions 

The Hodge decomposition theory plays an important role in the analysis of vector 
spaces in a 3D bounded smooth domain, our analysis on the boundary conditions will 
be based on this theory. To be self content, we give a simple L 2 version below. For 
more details, we refer |lll 14*0] and the references therein. 
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Let Q, C R 3 be a bounded smooth domain, H S (Q) denote the standard Hilbert space 
with H°(Q) = L 2 (Cl). Then the following estimate is well known. 



\\v\\ s < c(||V x v\\ s -i + ||V • v\\ s -i + \n ■ v\ s _i + |M| s -i) (2.1) 
for all v G H s (fl), s > 1 (see [2 [22]). 

Let u £ L 2 {9). Set 

u = v + Vi^ 9 
Note that V • u € // _1 (il). Let (p g solve 

Av? 9 = V • u in (2.2) 
ifg = on <9f) (2.3) 



It follows that 



Set 



Note that 



One has 



V • v = in n (2.4) 

DF = {ue L 2 (fl);V -u = 0} 
GG = {ue L 2 (n);u = Vcp,ip£ H%(tl)} 

(it, Vv?) =0, Vtt £ DF, y> £ #<}(fi) (2.5) 



Lemma 2.1 XTie following decomposition holds: 

L 2 (Q.) = DF © GG (2.6) 

Let u € -Di 7 . Then it • n is well-defined on dVt (see [23]) and 



Let if solve 

Set 
and 



u . n = / V • u = (2.7) 
an Jo, 



Aip = in n (2.8) 
5 n </3 = it • n on <9il (2-9) 



v = u — Vip 



H = {«£ L 2 (0); V • it = 0, in fi; it • n = 0; on d£l} 

DFG = {ue L 2 (Q);u = V<p, V • it = 0, / it • n = 0} 

Jon 

Hence v € H. Note that 

(«, v)=0, VueH,v £ DFG (2.10) 

It then follows that 

R 



Lemma 2.2 DF has the following decomposition: 

DF = H @ DFG 



(2.11) 



Note that u = V92 € DFG may not belong to the range of curl, and the range of curl 
V x H l {Q) is closed in L 2 (J7). DFG can be further decomposed to 

DFG = CG © HG (2.12) 

where 

CG = DFG n (V x H 1 ^)), HG = DFG n (V x H 1 ^)^ 

Let u = Vip G HG. Since 

= ((V x v),V(p) = / (nxw)-Vv9 (2.13) 
Jen 

for all u G i? 1 (r^), thus d T ip = 0, on with r being any tangential direction on c?0 
which implies (p is a constant on each component Tj of <9fi. So 

ilG = {V(/j; A99 = 0, <f = Ci on Tj} 

consists only smooth vectors, and is finite dimensional, which is called the harmonic 
gradient space. 

Remark 2.3 CG can also be expressed as 

CG = {u € L 2 (Q); u = V99, V • u = 0, / u ■ n = 0} 

Since CG C V x H (£1), we will call it curl type gradient space. 

Note that H DKer(Vx) is compact in L 2 (Q) due to (2.1). Set 

HH = HC\Ker{Vx) 

Then 

HH = {u G L 2 (n); V ■ u = 0, V x u = in fi 5 u ■ n = 0; on 917} 

This is called the harmonic knots space, which consists only smooth functions and is 
finite dimensional (see [LT]). Now H can be decomposed to 

H = FH © HH (2.14) 

where 

FH = H n (ifer(Vx)) 1 - 

In conclusion, we have 

Lemma 2.4 T/ie following decomposition holds: 

L 2 {£l) = FH®HH@CG@HG@GG (2.15) 



Then for any u £ L 2 (S1), it is uniquely written to 

u = Phhu + Pfhu + Pcgu + Phgu + Pggu (2.16) 
where Px denotes the projection on the corresponding subspace. 

It should be noticed that the space FH has the following expressions (see [9, 20, 36, 
47]). 

Lemma 2.5 The space FH can be expressed as 

FH = {ue L 2 (Q);V ■ u = 0, u ■ n = on dVt, F(u) = 0} (2.17) 

FH = {u; u = V x v, v E H\Q), V • v = 0, n x v = on dn} (2.18) 
where F(u) = means 

u • n = 

for any smooth cross section T,ofQ. 

It follow from (2.1) and the fact that HH C V x (FH n H 1 ^)) (see [TQ) respectively 
that 

Proposition 2.6 

L 2 (n) = V x [FH niJ^fi)) $FG©GG (2.19) 
Similarly, in general, it holds that 
Proposition 2.7 

F s (ft) = v x (fh n F s+1 (o)) e (hg n tf s (ft)) e (gg n H s (n)) (2.20) 

/or s > 0. 

It follows from ,$ZM ,$ZM and the fact that HH , HG are finite dimen- 

sional that 

Proposition 2.8 C°°(VL) n X is dense in H S (Q) Dl, s > for 

X = FH, HH, CG, HG, GG 

3 The Stokes operators 

In this section, we apply the Hodge decomposition theory to the Stokes problems with 
both the VSB and NSB conditions. We first consider a special Stokes problem with 
the VSB (3.1)-(3.3) and prove theorem 3.1.. Next, since the topology of the domain 
is assumed to be general, to avoid the uniqueness of the solutions for the general 
Stokes problems, we consider the perturbed Stokes problem associated with VSB ()3.25j) - 
(|3.27p . Based on theorem 3.1., by using the Hodge decomposition theory, we prove 
the associated Stokes operator is a self-adjoint extension of the associated positive 
definite bilinear form (see theorem 3.5.). The proof of theorem 3.5. is constructive, and 
the techniques can also be used to prove the well-posedness of the non-homogeneous 
problem (|3.55p - (|3.57|) (see theorem 3.7.). More generally, we can prove well-posedness 
of the boundary value problem (|3.63p - (|3.65p (see theorem 3.9.) by construction a 
contraction map. Finally, we identify the relationship between the NSB and VSB, and 
establish a similar theory for the Stokes problem associated with NSB (13.67p -( |3.69p . 
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3.1 A special Stokes problem 

Let us start by considering the following special Stokes problem with homogenous VSB 
condition 

-An = / in (3.1) 
V ■ u = in n (3.2) 
u ■ n = 0, nxVxii = 0on9n (3.3) 

with / G FH. Set 



W = {u G H 2 (n); n x (V x u) = on, 90} 



Then we have 



Theorem 3.1 The Stokes operator Ap = — A with the domain D{Ap) = W n F.ff is 
self-adjoint in the Hilbert space FH . 

Proof: It is clear that Ap = —A with the domain W n Fi? is symmetric. Since 
Cfi° (Q) (1 H is dense in H, it follows that j4p is densely defined due to the orthogonality 
of FH and HH and the compactness of HH. Let u G W. Since n x (V x u) = on 
<90, then -An = V x (V x u) G FiJ by lemma 2.5., thus A F maps W n Fff to Fif. 
Now, for any / G FH, it follows from lemma 2.5. that there is a $ G iif 1 ($1) satisfying 

V x $ = / in (3.4) 

V • $ = in n (3.5) 
$ x n = on dn (3.6) 

Due to proposition 2.7. and lemma 2.1., there is a v G FiJ n H 2 (Q) so that 

$ = V x v + Phg$ (3.7) 

Note that Phg^ x n = on dQ.. It follows that 

nx (V xv) = on dn (3.8) 

Then V x (P H g$) = and ([377]) imply that 

- At> = / in (3.9) 

Thus Ap ■ W n FIT — > F-ff is surjective. If / = 0, then integration by part shows 

||Vxu||=0 (3.10) 

It follows that u = due to the orthogonality of FH and HH and then Ap : WDFH — >■ 
Fii" is one to one. 

Noting that W and FH are closed in H 2 (£l) and L 2 (f2), and 

||Au|| < |H| 2 (3.11) 
we obtain from the Banach inverse operator theorem that 

Ha < c||Au|| (3.12) 
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The theorem was proved. 

Equivalently, we have shown the problem (3.1)-(3.3) has a unique solution u G H 2 (£l) 
for any / G FH. 

It follows from the proof of theorem 3.1. that 

Vx : H*(n) i y HF 

is also one to one and onto, where 

ifi(fi) = {u G fl^fi); V • u = 0; n x « = on dQ; (u,<p) = 0,V<p G HG} 

It follows from the trace theorem and the continuity of the divergence operator that 
is closed in H 1 ^). Then 

||u||i < c||V x u|| (3.13) 

follows from 

||V x u|| < (3.14) 
for any u G H^(Q). This yields immediately that 

Lemma 3.2 Let u G H^(Q). Then the following Poincare type inequality holds 

\\u\\ < c||V x u\\ (3.15) 

Let v G FH n ff^fi). Then there is u G ^(fi) such that V x u = v and 

(v,v) = (V x = (u, V x v) < \\u\\\\V x v\\ (3.16) 
This, together with (3.15), shows 

(v,v) < c||V x u||||V x n|| (3.17) 

Thus, one gets 

Lemma 3.3 Let u G FH n H 1 ^). Then the following Poincare type inequality 

\\u\\ < c||V x u|| (3.18) 

As a consequence, we can obtain 

Corollary 3.4 TTte operator Ap in theorem 3.1. is the self adjoint extension of the 
following bilinear form 

a(u, (f>) = (V x u, V x (/)), D(o) = V F = FH n fl" 1 ^) (3.19) 

in Ffl". 



Proof: From proposition 2.8., a(u,(j)) with D(a) = FH n iJ 1 (fi) is densely denned. 
Due to (2.1) and lemma 3.3., a(u,(f)) is closed and positive. It follows that there is a 
self-adjoint operator A with domain D{A) C D(a) such that 

a(u,<f>) = (Au,<f>),V<j> G Mflff^O) (3.20) 

for any u G D(A). It is clear that D(A F ) = W n Fff C D(vl) and Ait = -Au for any 
u G 1^(1 FH. Let u G and / = Au. It then follows that / G FiT. It follows from 

theorem 3.1. that there is a v G -D(^f) such that (3.1)-(3.3) are valid (with u replaced 
by v) and 

a(v,<f>) = (f,<i>) (3.21) 
for all 4> G Vf- On the other hand 

a(u,0) = (A«,0) = (/,</>) (3.22) 

for all cj) G Vf, hence 

a(«-»» = (Vx(u-^Vx^) = (3.23) 

for all 4> G Vf- Taking 4> = u — v shows that Vx(u-«) = 0. Thus u = v due to (2.14). 
Thus D(A) = D(A F ) and A = A F . 



Denote by V' F the dual space of V F respect to the L 2 inner product. Then the 
notation of weak solutions can be extended for / G V F : u is called a weak solution to 
(3.1)-(3.3) for / G V' F if 

a(u,<t>) = (f,cf>), V<f>eV F (3.24) 
3.2 The Stokes problem with VSB condition 

Next, we consider the Stokes problem with general VSB condition. Since the domain 
is allowed to have general topologe, the kernel of —A may be not empty. To avoid it, 
we consider the following boundary value problem instead: 

(J - A)u + Vp = / in n (3.25) 
V • u = in Q (3.26) 
u-n = 0, n x (V x u) = f3u on dtt (3.27) 

where f3 is a nonnegative smooth function. 
Define 

V = H^tyDH 
Wp = {u G H 2 (Q); n x (V x u) = flu on, dfl} 
Define a bilinear form as 

ap(u,4>) = (u,(f>) + ap(u,4>) 

where 

u,(f>)= I (3u-<t>+ I (Vxm)-(Vx^ (3.28) 
Jon Jci 



with the domain D(ap) = V. u G V is said to be a weak solution to the boundary value 
problem (13351) - ([33T1) on H for f € V if 

~ap{u,4>) = (/,<£), V0 in V (3.29) 
where V is the dual space of V. Based on theorem 3.1., we can prove 

Theorem 3.5 The self-adjoint extension of the bilinear form ap(u,4>) with the domain 
~D(ap) = V is the Stokes operator Ap = I + P(— A) with D(Ap) = Wp D H, and Ap is 
an isomorphism between D(Ag) and H with compact inverse on H . Consequently, the 
eigenvalues of the Stokes operator Ap can be listed as 

1 < 1 + Ai < 1 + A 2 > oo 

with the corresponding eigenvectors {ej} C Wp, i.e., 

Apej = (1 + Xj)ej (3.30) 

which form a complete orthogonal basis in H . Furthermore, it holds that 

(1 + AOIMI 2 <ap(u,u) < -l—\\Apu\\ 2 , VueD(^) (3.31) 

1 + Ai 

Proof: It is clear that a,p(u, <f>) with the domain D{ap) = V is a positive densely defined 
closed bilinear form. Let Ap be the self-adjoint extension of ap(u,(f>). It follows that 
Wp fl H C D(Ap) and Apu = u + P(— An), for any u G fl H by integrating by 
part. It remains to show that D(Ap) C fl H. Let it G P(^4/?) and / = Apu. Since 
D(Ap) C D(ajg) = I/, it follows from (3.29) that 

Mi<c||/|| (3.32) 

Let n(x) and f3(x) be internal smooth extensions of the normal vector f3 respectively. 
Then (5{x)u x n(x) G H l {Vi). Proposition 2.7. yields 

P{x)u x n(x) = V x v + V/i + Vg (3.33) 

with V/i = PhgW(x)u x n(x)), V# = P GG ((3(x)u x n(x)) and u G FH n # 2 (0). It 
follows that 

||Vs||i<cHi (3.34) 

since g satisfies 

Ag = V • (/3(x)u x n(x)) in fi (3.35) 

5 = on <90 (3.36) 

Since HG is finite dimensional, so 

\\P HG (P(x)u x n(x))||i < c||Phg09(z)« x n(x))|| < c||u|| (3.37) 

It then follows from (2.1) and lemma 3.3 that 

IMI2 < c||V x v\\i < c\\u\\i < c\\f\\ (3.38) 
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Integrating by part and noting that n x V/i = 0, n x Vg = on the boundary, we have 

/ (V x v) ■ (V x 0) + / fin x (it x n) ■ <t> = (-Aw, <f>) (3.39) 
Jn Jan 

for all G if 1 (ri). It follows from n x (u x n) = u on the boundary and the definition 
of the weak solution that 

(V x (u - v)) ■ (V x 4>) = (P FH (f -u + Av),4>) (3.40) 

for all <f> G H l {Vt) n FH". Note that Vxa = Vx ffcff (u) and Pfh(u) € H 1 ^) n FiT. 
It follows that 

a(P FH (M)-«,0) = (JVH(/-« + Ai;),^) J G H^tyHFH (3.41) 
It follows from theorem 3.1. that Pfh(u) — v G W and 

H-Pfh («) - «|| 2 < c(||/|| + ||Au|| + H|) (3.42) 
Since HH is finite dimensional, it holds that 

\\Phh(u)\\ 2 < c\\u\\ (3.43) 
One gets from (|532|l . (|g38|l . (|515|l and (13331) that 

N|2<c||/|| (3.44) 
Since Pfh{u) — v G W, it holds that 

nxVxu=nxVx Pfh(u) = nxVxv = f3u on dfJ (3.45) 
Thus we have shown u G H H . Integrating by part in (3.29) yields 

(u - Au - f, (f>) = (3.46) 

for all (j) G V, which implies 

u - Au + Vp = f in U (3.47) 

with p given by 

-Ap = in n (3.48) 
(Vp) • n = Au ■ n on dQ (3.49) 

It is noted that ||Agii|| is an equivalent norm of H 2 (Q) on Wr D i? due to (3.44) and 

HA^H < ||w|| + ||Au|| < c||it|| 2 (3.50) 

for all u G Wr D i?. The theorem was proved. 

Let V be the dual space of V with respect to the L 2 inner product, u G V is called 
a weak solution to (3.1)-(3.3) for / G V if 

a p (u,4>) = {f,<f>),y<l>€V (3.51) 

By using a standard density argument, one can show 
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Corollary 3.6 For any f £ V , the boundary value problem \3.25\) -i 3.21) has a unique 
weak solution u £ V 

Now, let b £ Hz{dQ) and b ■ n = on <9ft. From the extension theorem, it has an 
extension denoted by b(x) £ if 1 (ft). Similar to the proof of theorem 3.5., one can show 
that there exists a £ H 2 (Q) n FH such that 

n x (V x <J>) = 6 on <9ft 

It follows that $ solves the following problem: 

u - An + Vp = / in ft (3.52) 

V • u = in ft (3.53) 
n • n = 0, n x (V x u) = b on dft (3.54) 

with f = u + P(-A$) and Vp = A$ - P(A$). This fact and theorem 3.5. for /3 = 
yield 

Theorem 3.7 Zei 6 £ Hz(dQ), b ■ n = and A > 0. TTien £/te following problem 

An - An + Vp = / in (3.55) 

V • n = in ft (3.56) 
u ■ n = 0, n x (V x u) = b on 9ft (3.57) 

has a unique solution u £ H 2 (Q) for any f £ H. 

The boundary value problem (I3.55l) - (l3.57p also have a weak formulation 

\(u,<f))+ f (V xn)-(V X0)+ / b-<f> = (f,<f>), V<£inV (3.58) 

Jan 

Similar to corollary 3.4., one has 

Corollary 3.8 Lei 6 6 H~2(d£l), b ■ n = 0. T/ien /or any / £ 1/', the boundary value 
problem i3. 55\)-(3.57\) has a unique weak solution u £ V in the sense of i3. 58]) . 



We omit the details of the proof here, and refer to [26] for the definition of the weak 
tangential trace H~i(dQ). 

For any given smooth and nonnegative function /3, we define the map 

T : (ft) flfl'^Fc (ft) n H 

by u = Tv determined by (|3.58j) with b replaced by (3v+b and / = 0. Let V{ £ (ft)["W 
and Ui = Tvi, i = 1,2. It then follows from (3.58) that 

A||ni - n 2 || 2 + ||V x (m - n 2 )|| 2 + / /3(ni - n 2 ) • ( Vl - v 2 ) = (3.59) 

Note that 

/3(ni - u 2 ) ■ (ni - v 2 )\ < c\\ui - n 2 || i ||t>i - v 2 \\ i, , (3.60) 
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and 



It follows that 



Ml^i (n) < c IMIIMbi(tt) - c llvllll v x vll> V <P G v ( 3 - 61 ) 



\\ui-u 2 \\ 2 i < c\~h\vi -v 2 \\ 2 i (3.62) 

for A > 1. Take A large enough such that T becomes a contraction map on Hz (£2). It 
follows that 

Tv = v 

has a unique solution ^ on #2(0) and ^ = £ 

For any / € V*', let v be the weak solution of (^25> (13371) with / = / - (1 - A)f. It 
is clear that tt = v + 1 I / €yisa weak solution of the following problem: 

u - An + Vp = f in U (3.63) 
V • it = in n (3.64) 
u ■ n = 0, n x (V x u) = /3-u + 6 on dtt (3.65) 

in the sense that 

(u,4>) + / (V x u) • (V x 4>) + I (Pu + b)-<f>= (f,4>), V <f> in 7 (3.66) 

The uniqueness can be proved in the same way as for theorem 3.1.. We conclude 

Theorem 3.9 Let b € H~^(dQ,), b ■ n = 0. Then for any f € V , the boundary value 
problem 13. 63)) - {3. 65]) has a unique solution u € V in the sense of 13. 66]) . Moreover, if 
f £ H and b G H*(dQ), then u € H 2 (Q). 

3.3 The Stokes problem with the NSB condition 

We can establish a similar theory for the Stokes problem with the NSB just as with 
VSB. For completeness, we sketch it here. Consider the following Stokes problem with 
the NSB condition. 

(I - A)u + Vp = f in Q (3.67) 
V • u = in n (3.68) 
u ■ n = 0, 2(5(u)n) T = —^u T u on d£l (3.69) 

where 7 is a nonnegative smooth function. 
Define 

W-y = {«£ H 2 (n); 2{S{u)n) T = - 7 u T on, 3ft} 



and a bilinear form 

a 7 (u, 4>) = (u, <f>) + a 7 (u, </>), Z?(a 7 ) = V 



where 

a 7 («,0)= f 7^-0 + 2 / S(u)-S((f>) (3.70) 

and 5(u) • S(cfr) denotes the trace of the product of the two matrices. 

u is said to be a weak solution to the boundary value problem (I3.67p - (l3.69j) on H for 
/ € V if 

a>, $ = (/,<£), V<£ in V (3.71) 
where V is the dual space of V. 

To compare it with the VSB case, we first calculate that 
Lemma 3.10 Let u € H 2 {Q) and u ■ n = on the boundary. It holds that 

(2(S(u)n) -uxn) T = GD(u) T (3.72) 

with GD{u) = -2S(n)u. 
Proof: Note that 



and 



It follows that 



d n u = -u x n + S(u)n (3.73) 



d T u = ^ooxT + S(u)t (3.74) 



2(5(u)n) • r = <9 T -u • n + d n u ■ r (3.75) 
(n x w)r = <9 T ti • n — d n u ■ r (3.76) 



and 



2{S{u)n) ■ t + (n x w)t = 2d T u ■ n (3.77) 
Note that u ■ n = on the boundary. It follows that 

d T u ■ n = —u ■ d T n (3.78) 

We conclude that 



Note that 
thus 

Note that 
and 



(2S(u)n -ujxn)-T = -2u ■ d T n (3.79) 



9 T n = i(Vxn)xT + S(n)T (3.80) 



(2S(u)n - uj x n) ■ r = ((V x n) x u) • r - 2S(n)u ■ r (3.81) 

uxt = Xn (3.82) 

(V x n) ■ n = (3.83) 
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on the boundary. It follows that 

(25(«)n - OJ x n) ■ r = -2S(n)u ■ r (3.84) 

Set 

GD(u) = -2S{n)u 

The lemma is proved. 

It follows from a simple calculation and by using the density method that 
Lemma 3.11 Let u G H 1 ^) Pi H. Then 

2 [ S(u) ■ S((p) = [ (V x u) ■ (V x <f)) + [ GD(4>) ■ u (3.85) 
Jn Jci Jan 



GD{(f>) ■ u = [ GD(u) ■ <f> (3.86) 
Jan 



ian Jan 
As a counterpart of theorem 3.5., we can obtain 

Theorem 3.12 The self-adjoint extension of the bilinear form a 7 (u, <fi) with domain 
D(a 7 ) = V is the Stokes operator = I + P(— A) with D(Ay) = W" 7 PI H , and is 
an isomorphism between D(^4 7 ) and H with a compact inverse on H . Consequently, 
the eigenvalues of the Stokes operator A^ can be listed as 

l<l + Ai<l + A 2 > oo 

with the corresponding eigenvectors {ej} C W 7 , i.e., 

A 7 ej = (1 + Xj)ej (3.87) 

which form a complete orthogonal basis in H . Furthermore, it holds that 

(1 + Ai)||n|| 2 < ap(u,u) < - — \\A jU \\ 2 , Vu G D(A 7 ) (3.88) 

Proof: It suffices to show that D(A^) C W 7 H H since the rest is similar to the proof 
of theorem 3.5.. Let u G D(A 1 ) and / = A^u. Since D(A 1 ) C D(ap) = H l (£l) n H, it 
follows from (|3.7ip that 

Hi < c||/|| 2 (3.89) 

Let n{x) and 7(x) be internal smooth extensions of the normal vector n and 7. Then 
( / y(x)u + GD(u)) x n(x) G H 1 ^). Due to proposition 2.7., one has 

{j(x)u + GD(u)) x n(x) = V xv + Vh + Vg (3.90) 

with v G H 2 (Q) n F#, V/i = P HG ((j(x)u + GD(u)) x n(x)) and V# = P GG ((7(x)u + 
GD(u)) x n(x)). Similar to the proof of theorem 3.1., one can get 

1Mb < c|| V x < c||it||i (3.91) 
1fi 



Note that n x (V/i) = and n x (Vg) = 0. Thus 

/ (V x v) • (V x cp) + [ (<yu + GD(u)) ■ cj) = (-Av, G V (3.92) 

Jan 

Then the definition of the weak solution and lemma 3.11. imply 



(Vxu)-(Vx0) + / ju-<f> + GD(4>) ■ u = (f — u,(p), \A/> € F (3.93) 
n Jan Jen 

Combine them and note (|3.86p to get 

/ (Vx(ti- u)) • (V x 0) = (Pp H (/ - « + Av), 0), G 7 (3.94) 
Jn 

Note that V x u = V x P FH (u) and Pp(u) G iT 1 ^) n FH. It follows that 

a(P FH (u)-^0) = (P ra (/-« + A^,0) ! V0 G H^tyHFH (3.95) 
Since PfhU -u + Av) G FJT, so Pfh(u) -v£W, and 

||Pfj(w)-«|| 2 < c(||/|| + ||w||i) (3.96) 
Since -ff-ff is a finite dimensional, so 

||Phh(«)||2<c|H| (3.97) 
It follows from ([3T89]) . (^9Tj) . (^96]) and fllggD that 

|H| 2 <c||/|| (3.98) 

Note that 

(V x u) x n = (V x P F h{u)) xn = {Vxv)xn = -yu - GD{u) (3.99) 
It follows that 

2(S(u)n) T = ((V x u) x n + GD{u)) T = -ju T (3.100) 
The theorem was proved. 

Similar to the discussion for VSB, we have 

Theorem 3.13 Let b G H~2(dQ) ! b • n = 0, 7 be a nonnegative smooth function on 
the boundary. Then for any f G V , the following boundary value problem 

u- Au + Vp = / in n (3.101) 
V • u = in n (3.102) 
u ■ n = 0, 2{S{u)n) T = - 7 ti T + 6 on 90 (3.103) 

has a unique solution u G V in the sense that 

(u,<f>) + [ ( 7 u + 6) • </> + 2 [ S(u) • 5(0) = (/, <j)), in V (3.104) 
Jan Jn 

Moreover, if f G H and b G Hz(d£i), then u G H 2 (Q). 
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4 Functional setting of the LNS-a equation 



In this section, we formulate the following boundary value problem for the LNS-a 
system: 

d t v - Av + V x v x u + Vp = in tt (4.1) 

V • v = in n (4.2) 
u — aAu + Vp = v in f2 (4.3) 

V • u = in n (4.4) 

with the VBS conditions 

v ■ n = 0, n x V x v = (3v on <9f2 (4.5) 

« • n = 0, n x V x u = f3u on <9fi (4.6) 

Due to theorem 3.5., A a = I — aPA is also a positive definite self-adjoint operator with 
domain D(A a ) = Wp n H for any a > 0. We have 

Proposition 4.1 T/ie linear operator T a = A~ x : H i-> W 7 ^ n H is well defined with 
u = T a v G WJg n if given &?/ i/ie Stokes boundary value problem 

u - aAu + Vp = v in U (4.7) 

V • u = in Q (4.8) 

n ■ n = 0, n x (V x u) = fiu on dn (4.9) 

and is bounded, i.e. 

\\uh < c a \\v\\ (4-10) 
for some constant c a depending on a. 

We now estimate the nonlinearity. Let v G V C so that T a t; is defined. Set 

B(v,u) = P(V xv xu), Vu€Wp, veV (4.11) 

B a (v) = B(v,T a v), veV (4.12) 

for a > 0. Then we have 

Lemma 4.2 TTte nonlinearity B a (v) : V ^ H is locally Lipshitz for a > 0. 

Proof: Clearly, B a is well-defined due to (4.10). For any v±, v 2 G V, 

||5 a («i) - B a («2)|| < || V x (m - «2) x T Q vi - V x v 2 x T a ( Vl - v 2 )\\ (4.13) 
Note that for all (f> G V, V G 

||Vx(0)x^||<c|H|i|H| L =o (n) (4.14) 

and 

< c||io||i||ii;||2 (4-15) 
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It follows that 

\\B a ( Vl ) - B a (v 2 )\\ < cdl^Hx + ||«2||i)||(«i - u 2 )||i (4.16) 
which implies the lemma. 

We now can formulate the initial boundary problem of the LNS-a equations (4.1)- 
(4.6) as an abstract equation 

v' - PAv + B(v,u) = (4.17) 
u = T a v (4.18) 

on H, with a parameter a G (0, oo). 

The weak solutions of the initial boundary problem can be defined as below. 

Definition 4.3 (v,u) is a weak solution of (4-l)-(4-6) with a > for LNS-a equations 
with initial data vq G H on the time interval [0,T) if v G L 2 (0,T;V) n C w ([0,T]; H), 
v' G L 1 (0,T;V / ), u G L 2 (0, T; F n i? 3 (Sl)) n C w ([0, T\;Wp), v! G £*((), T;F) suc/i f/iai 

(t/, iu) + ap(y, w) + u),w) = 0, a.e. t G [0, T) (4.19) 
u = T a v, a.e. t G [0,T) (4.20) 

for all w £ V. 

For the special case a = 0, we define also the corresponding weak solutions for the NS 
as follows 

Definition 4.4 (v,u) is a weak solution of (4-l)-(4-6) with a = (NS equations) 
initial data vo £ H on the time interval [0, T) ifv,u G L 2 (0,T;V) n C w ([0, T]; H) and 
v',v' G L l {0,T;V) such that 

(v', w) + ap(v, w) + (B(v, u),w) = 0, a.e. t G [0, T) (4.21) 
ii = t;, a.e. t G [0,T) (4.22) 

for all w G V. 

For later use, one can also define the fractional powers of the operator Ap = I — PA 
in theorem 3.5., : D(A S ^) i-> H for s > by 

oo 

^« = E(l + A,-) a «je j (4.23) 
i=i 

for u G X^fci n j e i G -D(A|), where 



oo 



D(A^) = {u = Y J U j ey, ^(1 + A,) 2 * |u/ < oo [> (4.24) 



equipped with the graph norm 

■'lD(A-) 



= (A%u,Alu). 
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It can be checked easily that A s p : D(A s ^ rt ) i-> D(A t /3 ) is an isomorphism for all s, t > 0, 

i 

D(A 1 j3 ) = D(Ap) = ff 2 (n)n^, and D(A£) = V with equivalent norms \\u\\ 1 and 

) 

.ff 1 (fi)-norm. Denote by D{A~j 3 s ) the dual space of D(Ap) for any s > 0. Then the 
operator A^ can be extended to an operator: H i-> D( J 4^" S ) by 

(A£ u, u) = (u, ^ u), Vueff, w € D(A^) (4.25) 

It follows from the definition that 

iMi 2 D(A r) = iiv^ii 2 ' yu g ( 4 - 26 ) 

and A s : D(A^ +t ) i-> D(^4Jg) is an isomorphism for s, t G E, and furthermore, 

\\A^u\\ 2 = {A s pU,A p u) < \\A s pu\\ \\Apu\l Vu G D{A s p ) C\D{A p ) (4.27) 
holds true for all s, t G E. 

5 Well-Posedness of the LNS-a Equations 

In this section, we investigate the well-posedness of the initial boundary value problem 
of the LNS-a equations (4.1)-(4.6) by a Gelerkin approximation based on the orthogonal 
basis given in theorem 3.5.. 

5.1 Local well-posedness 

We start with the following local well-posedness result. 

Theorem 5.1 Let vq G H and a > 0. Then there is a time T* = T*(vq) > such that 
the problem (4-l)-(4-8) has a unique weak solution of (v,u) with initial data v$ on the 
interval [0, T*) in the sense of definition 4-1 for any T G (0,T*), which satisfies the 
energy equation 

^\\v\\ 2 + 2ap{v,v) + (B(v,u),v) = 0, on[0,T] (5.1) 
in the sense of distribution. Furthermore, if vo G V, then 

v G L 2 (0, T; Wp H H) D C([0, T);V) (5.2) 
v' G L 2 (0, T; H) (5.3) 



and the energy equation 

d_ 

dt 

is valid. 



d ap(v, v) + 2||PAv|| 2 + 2(B(v), -Av) = (5.4) 
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Proof: Let vq € H. Consider the following system of ordinary differential equations 

v' j (t) + X j v j (t)+g j (V) = (5.5) 
Vj(0) = (u ,ej) (5.6) 

j = 1, • • • m, where V = (vj) and 

g i (V) = (B(2?v j e j ,u m ),e i ) (5.7) 
u m = T a {Y.?Vjej) (5.8) 

Note that all norms are equivalent in a finite dimensional linear space. It follows from 
lemma 4.2. that (gj(V)) is locally Lipshitz in V and thus the systems is locally well 
posed and equivalent to the following partial differential equations 

v' m (t, x) - PAv m {t, x) + P m B(v m , u m ){t, x) = (5.9) 

u m = T a {Y.fvjej) (5.10) 

v m (0) = P m (vo) (5.11) 

where v rn (t,x) = T,™Vj(t)ej(x), and P rn is the orthogonal projection of H onto the 
space spin{ej}™. 

Taking the inner product of (5.9) with v m and noting that 

(P m B(v m ,u m ),v m ) = / V x v m x u m ■ (v m )dx (5.12) 
Jq 

one can get 

— \\v m \\ 2 + 2ap(v m ,v m ) + (B(v m ,u m ),v m ) = (5.13) 
It follows from the definition of T a that 

\{B(v m ,u m ),v m )\ < c\\v m \\i\\v m \\\\u m \\L°°(n) < c||v m ||i||f m || 2 (5.14) 

Note that 

U\\ 2 < cU\\l < c(Uf + a^A)) (5-15) 
for all (j) G V. It follows that 

^\\ v m\\ 2 + ap(v m ,v m ) < c(||w m || 2 + l)||w m || 2 (5.16) 
Hence, there is a time T > such that 

{v m } is bounded in L°°(0,T; H) 
{v m } is bounded in L 2 (0,T; V) 

Note that for 4> £ V, 

\(A p v m ,<l>)\ < \{v m ,(/))\ + \a p (v m ,<f>)\ (5.17) 

which implies that 

{Apvm} is bounded in L 2 (0, T; V') (5.18) 

Since 

||«m||L°°(n) = c\\T a v m \\ 2 < c\\v m \\ (5.19) 
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it follows that 

\{P m B(v m ,u m ),4>)\ = |(V x v m x u m ,P m 4>)\ < C\\v m \\i\\v m \\\\<f)\\i (5.20) 
for all (f) E V, which implies that 

{P m B(v m ,u m )} is bounded in L 2 (0, T; V') (5.21) 

Hence 

{v' m } is bounded in L 2 (0, T; V') (5.22) 

By using a similar argument in |20| . it shows that there is a subsequence also denoted 
by v m and a v G L°°(0, T; #) n L 2 (0, T; V) such that 

v m ^v in L°°(0, T; fl) weoA: - star (5.23) 
j) m — > v in L 2 (0,T;V) weakly (5.24) 
Vm — > v in L 2 (0,T;H) strongly (5.25) 

Consequently, u m = T a v m has the property: 

u m -> T a v in L°°(0, T; n JET) wea/c - star (5.26) 
u m -> T a t; in L 2 (0, T; 7 n H 3 (tt)) weakly (5.27) 
n m -> T a v in L 2 (0, T; n JET) strongly (5.28) 

Passing to the limit of a subsequence, it is showed that (v, u) is a weak solution in the 
sense of definition 4.3.. It also follows that the energy equation 

j t \\v\\ 2 + 2ap{v,v) + (B a (v),v) =0 (5.29) 
is valid on the interval [0, T] in the sense of distribution. 

Let v\ and V2 be any two solutions. Then w = V\ — V2 satisfies the following equation 

w' - PAw + P(B a ( Vl ) - B a (v 2 )) = (5.30) 
u>(0) = (5.31) 

and the energy equation 

iMI 2 + 2ap(w, w) + (B a ( Vl ) - B a (v 2 ),w) = (5.32) 
dt 

It follows from the local Lipshitz continuity stated in lemma 4.2. and the Gronwall 
inquality that 

\\wf < c(T)\\w(0)\\ 2 on[0,T] (5.33) 

which implies the uniqueness of the solution. Consequently, the convergence of the 
whole sequence follows. 

By the standard continuation method, there is a T* > such that the weak solution 
does exist on [0, T] for all T < T*, and if T* < oo then 

\\v(t)\\ -» oo, as t -4 T* 



9.0. 



Let Vq G V. Taking the inner product of (5.9) with — PAv m and noting that 

(PmB{v ), -PAv m ) = (B(v m , um), -P&v m ) (5.34) 

one gets 

dt 



d -Mv m ,v m ) + 2\\ - PAv m \\ 2 + (B a (v m ),-PAv m ) = (5.35) 



It follows that 

jj.ap(Vm,Vm) + || - PAv m \\ 2 < ||Wm||i||«m||ioo(j2) (5.36) 

Due to 

\\Um\\L°°(n) < Ca\\v m \\ (5.37) 

vq € V, the bounds of v m in L 2 (0,T;1/), and the GronwalPs inequality, one has 

{v m } is bounded in L°°(0, T; V) 

{v m } is bounded in L 2 (0, T; Wp D H) 

which, together with the uniqueness, implies that the whole sequence indeed converges 
in the sense 

v m ^fvin L°°(0, T; V) weak - star (5.38) 
v m -> v in L 2 (0, T; n H) weakly (5.39) 
i> m ^ w in L 2 (0, T; V) strongly (5.40) 

This completes the proof of theorem 5.1.. 
5.2 Global well-posedness 

Now, we prove the following global well-posedness result. 

Theorem 5.2 Ifvo£H,a> 0, then the solution v obtained in theorem 5.1. is global, 
i.e., T* = T*(v ) = oo. 

Proof: Let v be the weak solution on the interval [0,T]. Then, 

u = T a v € L°°(0,T\WpnH) (5.41) 
Taking u as a test function yields 

(«', u) + o^u, u) + u), u) = (5.42) 
Since v = {I — aPA)u, then 

2(t/,u) = ^(Huf + cra^,^) (5.43) 

in the sense of distribution on [0, T\. Note that 

(B(v, u),u) = / (V x v) x it • u = (5.44) 
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It follows that 
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+ aap{u, u)) + 2( / (3u ■ v + / (V x v) ■ (V x u)) = (5.45) 
Jan is) 

Due to the smoothness and the boundary condition for u, it holds that 

/ (V x v) ■ (V x u) = - I pu-v+ (-Ait) • v (5.46) 



Consequently 

It follows that 
and 



^(||u|| 2 + aap(u,u)) + 2(a /3 (u,u) + q||PAu|| 2 ) = (5.47) 



(|H| 2 + aap(u,u)) < (\\u \\ 2 + aap(uo,Uo)) (5.48) 

/ (ap(u,u) + a\\PAu\\ 2 )dT < (\\u \\ 2 + a^(u , «o)) (5-49) 
•/ o 

On the other hand, it follows from the energy equation (5.1) and a similar argument 
as for (5.17) that 

^IMI 2 + ap{v,v) < c\\v\\ 4 + 1 (5.50) 

Noting that 

IMI 2 < c(||u|| 2 + a 2 ||PAn|| 2 ) (5.51) 
it follows that ^ 

\\v\\ 2 + / ap(v,v) < c (5.52) 
Jo 

for some constant c depending only on vq and a. Thus T* = oo. The theorem is proved. 



6 Vanishing a Limit and the NS Equations 

In this section, we investigate the vanishing a limit of the solutions of the LNS-a 
equations (a — > 0) to that of the NS equations. We will prove both weak and strong 
convergence results. Then, the global existence of weak solutions and the local unique 
strong solution to the NS equations with the VSB condition are followed. 

6.1 Weak Convergence and Global Weak Solutions of the NS 

We first prove 

Theorem 6.1 Let vq € H, and (v a ,u a ) be the global weak solution stated in theorem 
5.2. corresponding to the parameter a > 0. Then for any given T > there is a 
subsequence u a i of u a and a (v°,u°) satisfying 

v° eL 2 (0,T;V)DC w ([0,T];H) (6.1) 
(v )' €L*(0,T;V') (6.2) 
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such that 

v a J v° in L 2 (0, T; H) weakly (6.3) 

V 



v a i -)> u° in L 2 (0, T; D(Ag 4 )) strongly (6.4) 



v° in L 2 (0, T; Vp) weakly (6.5) 

u 01 * -> v° in L 2 (0,T;D(A~*)) strongly (6.6) 

Moreover (y°,v°) is a weak solution of the initial boundary problem of the NS equations 
(4-l)-(4-6) with a = and satisfies the energy inequality 

^-\\v \\ 2 + 2ap(v ,v°)<0 (6.7) 

Proof: Let vq G H, T > 0, and (v a ,u a ) be the global weak solution to (4.1)-(4.6) 
corresponding to 1 > a > 0. It follows from ([5.47P that 

\\u a f + aa p {u a ,u a )+ [ {ap(u a ,u a ) + a\\PAu a \\ 2 )dT < c (6.8) 

Jo 

for some constant c independent of a. For any 4> € Wp n H, we have 

(B(v a , u a ), <j>) = [ (V x v a x u a )0cfa = 1 + 11 (6.9) 



where 

7 (n x v a ) ■ (u a x cp)dS 



dCl 

11= [ v a ■ (-u a ■ V0 - <j> ■ Vu a )dx 
Jn 

Since u ■ n = and (j) ■ n = on the boundary so 

u a x 4> = An on <9f2 

Hence 

/ = 

To estimate II, we note that 

V Q • (u a ■ VcP)dx\ < c(\\u a \\ + a\\PAu a \\)\\u a \\ L s {a) \\V<l>\\ L e {u) 

2 



W a \\l m) < c\\u a Wm < c\\u a \\l(\\u a \\ + ||PAu a ||)i 



\\\7<P\\ L e m < c\\Ap<j>\\ 

Then, due to (6.8), it holds that 

| / v a ■ (u a -V(f))dx\ < c((ap(u a ,u a ))^ + a\\PAu a \\ + a\\PAu a \fi)\\Ap(f)^ 
Next, 

/■(f Vu a )dx\ < c{\\u a \\ + a||PA« a ||)||u a || 1 ||0|| £ «» ( n) 



6.10) 
6.11) 

6.12) 
6.13) 

6.14) 

6.15) 
6.16) 

6.17) 
6.18) 



which implies that 

| / v a -((l)-Vu a )dx\<c((ap(u a ,u a ))^ +a\\PAu a \\+a\\PAu a \\l)\\A /3 <P\\ (6.19) 
Then for a < 1, 

\(B(v a ,u a ),(f>)\ < c(l + (a/3{u a ,u a ))^ +ai\\PAu a \\%)\\A l3 (j)\\ (6.20) 

It follows from (j6.8j) and ()6.20p that B(v a ,u a ) and then ^{v a ) are uniformly bounded 
in Ll(0,T;D(A^)). It follows from (4.7)-(4.9) that 

(l-a)v? + aA f) (v?)=v?, 

which yields immediately 

(l-a)\\A/u?\\ 2 + a\\A/u?\\ 2 = \\A/v?\\ 2 . 

Then 

||^ 1 <|| 2 <2||A^<|| 2 

for < a < This shows that tt Q are uniformly bounded in Lz (0, T; D(Aa )) as 
9t!) a are. Note that (|6.8p also implies that (u a ) are uniformly bounded in L 2 (0,T;V) 

and the duality between V = D(Al) and D{A~^ 1 ) with respect to the inner product of 
D(A~*), i.e., 

(A-*u,A-*<t>) = (Alu,Ap 1 <l>) 

By using the standard compactness argument (see [2H [20] ) , one can show that there 
exist a subsequence u a -> of u a and a v° such that 

n Q J v° in L 2 (0, T; V) weakly (6.21) 
u a i v° in L 2 (0,T;D(A^)) strongly (6.22) 



Note that 

\(B(v a ,u a ) - B(v°,v°),(j))\ < I + 11 (6.23) 

where 

/ = \(B(u a -v°,u a ) + B{v°,u a -v°),<P)\ (6.24) 

II = a\(B(PAu a ,u a ),</))\ (6.25) 
Similar to (j6.9[) and ()6.13j) . integrating by part yields 

\{B{u a — v°,u a ),(f))\ = \ [ (u a -v°)-{u a -V<j) + 4>-Vu a )\ (6.26) 

Jn 

Note that 

r (u a - v°) ■ (u a ■ V0)| < c\\u a - ^ O ||||n Q ||i||n Q ||J 6(n) ||V0|| L6(n) (6.27) 



and 

.a „ 1 1 2 ^ ii a „.0i it a 



u a -v"\\' <c\\u a -v u \\ _i\\u a -v u \\ i (6.28) 



9P, 



\\u a -v°\\ 2 ! < c\\u a -v°\\\\u a -t; ||i 
This, together with (j6.8|) . shows that 

(n Q -t; o )|-(u a -V0)<c||u a -t; o p _ x ||u Q - u || J \\u a \\ \ \\A p <j>\\ 



Hence 



(u a - v°)\ ■ {u a • V0) < c\\u a - u°p _ x (||«°||J + ||t; ||J)p^|| 



While 

( ," ,,0 , ( , . V" ,,'M i ^ il ,," ,,0 I ii ,," 

It follows that 



( u «- w °).(^.Vu Q )| < ||« a -u o ||||« a ||i||0||L«(n) 



(u a -v u )\- (u a ■ V<f>) < c\\u a - v u \\ 2 (||u a ||< + 11^11^)11^11 



Then 



\(B(u a — V ,U a ),(f>)\ < c\\u a — B f _j (l+||u a ||J + ||«lJ)||480|| 

Similarly, one can obtain 



\(B(v°,u a -v°),(f>)\ <c\\u a -v°\\* i (l+K||f + ||i>°||?)ll^ll 



It follows that 



I < C\\u a - V°\\' ! (1 + IKIlf + ||^ ||f )\\Ag(t>\\ 
D(A*) 



*i8 

Similarly, 

|(5(PAn Q ,n Q ),0)[ = [ / (PAu a )-(u a ■ V 4> + cf) ■ V u a )\ 

Jn 

Then 

\(B(PAu a ,u a ),<t>)\ < c||PA-u a ||||u a ||i||^0|| 

It follows that 

// < ca^(a\\PAu a \\ 2 + |K|| 2 )||</>|| 2 
It follows from (16^31) . (|Q51> . (Q9l) . JEHU and ([OTP that 

B(v a i,u aj ) -> 5(v ,u ) m L^CT;!)^ 1 )) strongly 

which enables us to pass the limit in (4.19)-(4.20) to show that v° satisfies 

((v )', <P) + ai3 (v°, <j)) + ((V x v°) x v°, </>) = 0, a.e. t 

for all (p G C°°(U) CiV in the sense of distribution on [0,T]. Note that v° G L 2 
implies (v )' G £3(0,T; V'). Thus (IQTj) is also valid for all (f> G F. 
Due to (pTTTj) . it holds that 

-^-(IKH 2 + aa p {u a , u a )) + 2ap{u a , u a ) < 
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Passing to the limit and noting the weak lower semi-continuity of the norm, one gets 

j t \\v \\ 2 + 2ap(v ,v°)<0 (6.43) 

Note that 

(„« - v°, 4>) = {u a - v°, <j>) + a((A*u a ,A}<t>) - (u a ,<p)) (6.44) 



i 



for <j> € -D(A /3 4 ). Then 



\v a -v°\\ 2 _! <|K-?; || 2 _i +ca5( a ||PAM Q || 2 + |K|| 2 ) 



It follows that 



v a i ->• w° in L^T;!}^*)) strongly (6.45) 



Note that 

(w Q - t; , 0) = (u a - v°, 4>) - a(PAu a ,(f>) (6.46) 
It follows from ([621) and <t6^TTl that 

t; * -»■ v° in L 2 (0, T; JT) weakly (6.47) 
Hence, the theorem is proved. 

6.2 Strong Convergence and the Strong Solutions of the NS 

We now turn to the strong convergence of the strong solutions of the LNS-q to that of 
the NS equations, and prove 

Theorem 6.2 Let vq € V and (v a ,u a ) be the strong solution stated in theorem 5.1. 
corresponding to the parameter a > 0. Then there is aT > and a v° in L°°(0, T; V)f\ 
L 2 (0, T; Wp n H) such that 

v a -> v° in L 2 (0, T; Wp n H) weakly (6.48) 
v a -> v° in L 2 (0, T; V) strongly (6.49) 

u a -> u° in L 2 (0, T; % n JET) weakly (6.50) 
-u a -> v° in L 2 (0, T; V) strongly (6.51) 

loif/i t>° freing a weak solution to the initial boundary problem of the NS equation (4-1)- 
(4-6) with a = which is unique and thus called the strong solution. Consequently, it 
can be extended to the maximal existence time interval [0, T*) such that ifT* < oo then 

\\v ||i — > oo, as t — > T* 

Moreover, the following energy equation holds: 

^ap(v°, v°) + 2||PAt; || 2 - 2(B(v°, v°), PAv°) = (6.52) 
at 



Proof: It follows from the energy equation (5.5) that 



j t a p {v a ,v a ) + HPA^I 2 < c\\B(v a ,u a )\\ 2 (6.53) 



Note that 



\\B(v a ,u a )\\ 2 < c j |V x v a \ 2 \u a \ 2 dx < c||V x v a \\ 2 L3in) \\u a \\ 2 L6(n) (6.54) 

||V x v a \\ 2 L3{n) < c(\\v a \\ + UPA^IDII^Il! (6.55) 

\\u a he(n) < c||ti a ||i (6.56) 

and 

IK||i < c||v a ||i (6.57) 

which follows from the fact that 

\\u a \\ 2 + aap(u a ,u a ) = (v a ,u a ) (6.58) 

ap(u a ,u a ) + a\\PAu a \\ 2 = ap(u a ,v a ) (6.59) 

Consequently, 

^-ap(v a ,v a ) + -\\PAv a \\ 2 < c(l + |K|| 2 )IKII? (6.60) 
Combining this with similar estimates for (5.17) yields 

|(K|| 2 + ap(^,v a )) + \(\\v a \\ 2 + llPA^f) < c(l + ~ap(v a ,v a )) 3 (6.61) 
Comparing it with the following ordinary differential equation 

with y(0) = ap(vo,vo) shows that there is a time T such that 

v a is uniform bounded in L°°(0, T; V) n L 2 (0, T; Wp n H) 

It follows from this, (6.55)-(6.60), and (4.19) that 

(v a )' is uniform bounded in L 2 (0,T; H). 

Hence, by using the standard compactness argument, we find a subsequence v a i of v a 
and a v° such that 

v a * -»• v° in L 2 (0, T; Wp n if) (6.62) 

t/^ ->• w° m L 2 (0,T; V) strongly (6.63) 

which enables one to pass to the limit to find v° € C([0, T]; V) PI L 2 (0, T; n iJ) such 
that (v°,f°)is a (strong) solution of the NS equations. 

Let v\ and v® be two strong solutions to the Navier-Stokes equations with same initial 
data. Set w = v® — v!>. Then 

^H| 2 + 2ap(w, w) + 2(B(t>°, v °) - B (tM), «,) = (6.64) 



Note that 

\(B(vl n?) - B(v° 2 , v° 2 ),w) | < \(B(w, u?),u>)| + | (B(v° 2 ,w),w) \ 

< a p (w,w)+c(\\v^\\ Loo (t) + \\V xn°|| 4 )|H| 2 

which, together with (5.53) and GronwalFs inequality, yields ||to|| = 0. Thus we have 
obtained the uniqueness of the strong solution to the initial boundary value problem for 
the Navier-Stokes equations. By the standard continuation method, the strong solution 
can be extended to the maximum existent time interval [0, T*) D [0, T], and the energy 
equation follows from the smoothness of the solution. Consequently, the convergence 
of the whole sequence of v a follows. 



Finally, we prove the convergence of u a . It follows from (j4.3j) that 



V x u a - aA(V x/) = Vx v a , in n 



(6.65) 



Taking the inner product of above equality with — A(V x u a ) and integrating by part, 
we can get 



Au a \\ 2 + a||(Vxfn Q || = (Au a ,Av 



3„.ai 



Jan 



v a - u a ) 



(6.66) 



To handle the last term on the right hand side above, we use the fact v a — u a 
n x ((i> Q — u a ) x n) on dtt and the Stokes formula to get 



<><> 



Au a ■ ((3(v a - u a )) 



Au a ■ (n x {P{v a - u a ) x n)) 



<)<> 



(n x An") • (/3(n Q - v a ) x n) 



/ 



(V x (Au a )) ■ ((3(u a -v a )xn)- / Au a ■ V x ((3(u a - v a ) x n) 



where we have extended f3 and n smoothly to Vt. It follows from (4.7) that 
\\ v a - u a \\ 2 = {-aAu a ,v a -u a ) < a\\Au a \\ \\v a - u a \\ 

which yields 



,.o „,a 



u a \\ < a\\Au a \ 



It follows from (6.57) and (6.68) that 

An Q • V x ((3(u a - v a ) x n) 



< -||An a || 2 + c|K||? 



for suitably small a. Using (6.68) again gives 

|/ n (V x Au a ) ■ (f3(u a - v a ) x n)\ 

< -a [ |V x {Au a )\ 2 + a- 1 c\\u Q -v a \\ 2 
2 Jn 



< -a [ \V x (Au a )\ 2 + ca\\Au a \\ 2 
2 Jn 



Collecting (6.66), (6.67), (6.69) and (6.70) leads to 

||An a || 2 + a||(Vx) 3 n a || 2 < c\\v a \\l 



(6.67) 



(6.68) 
(6.69) 



(6.70) 



(6.71) 



for suitably small a. This, together with the bound of dt u a in H, implies the desired 
convergence in (6. 50), (6. 51). Thus the theorem is proved. 
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6.3 Estimates on Convergence Rates 



Finally, we study the rates of convergence in the case of strong solutions. We start with 
the case that the limiting Navier-Stokes system has a strong solution. 

Theorem 6.3 Let vq G V and v° be the strong solution to the Navier-Stokes equation 
with initial data vq on any given finite interval [0, T] with T > 0. Then there exists a 
cto > such that for each a G (0, ceo], the LNS-a with the initial data vo has a unique 
strong solution (v a ,u a ) on the same interval [0, T] satisfying 

sup \\(v a ,u a )-(v°,v°)\\ 2 + [ T \\(v a ,u a )-(v°,v°)\\l(t)dt<ca (6.72) 

0<t<T JO 

sup \\v a -v°\\ 2 + [ \\v a - v°\\ 2 2 dt < cah (6.73) 

0<t<T Jo 
with c being a positive constant depending on v° . 

Proof: Thanks to the local well-posedness of the strong solution to the initial-boundary 
value problem for the LNS-a and the standard continuation arguments, theorem 6.3. 
will follow immediately from the following a priori estimates. 

Proposition 6.4 Let T\ G (0, T] and (u a ,v a ) be the strong solution to the LANS-a 
with the initial data v° on the interval [0, Ti] with the property that 

\\v a \\j(t)+ f \\v a \\ 2 2 {T)dT < c foriG[0,Ti] (6.74) 
Jo 

with a positive constant cq depending only on T and v° G L°°(Q, T : V) H T T (0, T : Wp). 
Then there exist uniform constants a\ and c with the same dependence as Co such that 

sup ||(^,n«)-( W °,t; )|| 2 + [ TL \\(v a ,u a )-(v ,v )\\l(t)dt<ca (6.75) 

0<t<7i JO 

sup \\(v a - v°)\\j(t) + \\v a - v°\\ 2 2(t)dt < ca 1 * (6.76) 

0<t<T! JO 

for a G (0, ai]. 

Remark 6.5 Assuming Proposition 6.4- for a moment, one can verify the a priori 
assumption (6.74) by choosing 

f T 

c = 1 + 4 sup ||v ||f + / \\v°\\ldt 

0<t<T Jo 

and using (6.75) and (6.76) to choose a®. Thus (6. 75), (6. 76) hold for T\ = T. This 
yields theorem 6.3. immediately. It remains to verify proposition 6.4- 
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Proof of Proposition 6.4. Set w = v a - v°. Then it holds that for t G [0,Ti], 

^\\w\\ 2 (t) + 2ap(w, w)(t) + 2(B(v a , u a ) - B(v°, u°),w)(t) = (6.77) 

Note that 

2(B(v a , u a ) - B(v°, v°), w) = 2(B(v a ,v a ) - B(v°, v°),w) + 2a((V x v a ) x P(Au a ),w) 

\(B(v a ,v a ) - B(v°,v°),w)\ = \(B(w,v°),w)\ < c||v || 2 |Hli HI 
|(V x v a ) x P(Au a ),w\ < c\\P(Au a )\\ \\v a \\ 2 |H|i 
It follows that for all t G [0, T\] 

j t \H\ 2 (t) + ap(w,w)(t) < c(l + \\v°\\ 2 2 )\\w\\ 2 + ca 2 \\P(Au a )\\ 2 \\v a \\ 2 (6.78) 

Due to (6.59), (6.60) and (6.74), one has 

a||P(Au Q )|| 2 < cc for all t £ [0,71] (6.79) 
which, together with (6.78), yields 



— \ 



w\\ 2 {t) + ap(w,w)(t) < c(l + || ||^|| 2 (t) + ccoa||v a ||!(t) (6.80) 



Since w(0) = 0, so Gronwalhs inequality leads to 

sup < t < Tl IHi)|| 2 + /o Tl IHIiW* 

< (cacoe^o^a+lKHK*))^) /^IKHK^t (6.81) 

■/ o 

< ca c\e CCx = C2« 

with ci depending only T and the L 2 (0,T : 77 2 )-norm of v°. 

To prove (6.75) for u a , we note that (4.3) implies that 

\\u a -v°\\ 2 + aafs{u a ,u a - v°) = (v a - v°, u a - v°) (6.82) 

ap{u a - v°,u a - v°) + a(P(Au a ),A(u a - v )) = a p {v a - v°, u a - v°) (6.83) 
It follows from (6.82), (6.83), and (6.74) that 



\u a -v°\\ 2 + a/3 (u a -v°,u a -v°) < \\v a -v°\\ 2 + a ai3 (v°,v ) 

< C2CX + CCoCt 



Since 



(6.84) 



a{P{Au a ),A{u a -u )) = a\\P(Au a )\\ 2 - a(P(Au a ), Av°) 

> i«||P(An-)|| 2 -|||A,; || 2 

This, together with (6.84), shows that 

a/3 {u a - v°,u a - v°) + a\\P(Au a )\\ 2 < ap(v a - v°,v a - v°) + a||A-u°|| 2 (6.85) 
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Hence, one obtains from (6.82), (6.84), and (6.85) that 

sup \\u a -v°\\ 2 + / \\u a - v°\\\dt < c 3 a 

0<t<Tx Jo 

This completes the verification of (6.75). It remains to prove (6.76). 

By the definition of strong solutions, one has that for a.e. t G [0,Ti], 

jap{w, w) + 2\\PAw\\ 2 + 2(B(v a , u a ) - B(v°, v°), PAw) = (6.86) 

Rewrite the last term on the left hand side above as 

2(B(v a , u a ) - B(v°, v°), -P(Aw)) = 1 + 11 
One can estimate each term as 

= \2(B(v a ,v a ) - B(v°,v°), -P(Aw))\ < c(\\v a \\ 2 + |H| 2 )|Hli \\ PA ™\\ 

\II 2 \ = |2a((V x v a ) x (PAu a ),PAw)\ < c^PAu^ \\v a \\\ \\v a \\\ \\PAw\\ 
It follows that for a.e. t G [0,71], 

j t ap{w,w){t) + \\PAw\\ 2 {t) 

< c(l + \\v°\\l + \\v a \\ 2 2 )ap(w,w) ( 6 - 87 ) 
+ca 2 \\PAu a \\ 2 ||« a ||i \\v a \\ 2 + ca 2 \\PAu a \\ 2 \\v a \\ 2 

As a consequence of (6.57), (6.74) and 

aWPiAu )^ = IK-t^Hi 

one has 

ca 2 ||PAu a ||i ||t; Q ||i \\v a \\ 2 

< c\\v a -u a \\\ ||v a ||i \\v a \\ 2 (6.88) 

< c||f a ||i |H|f ap(w,w) + c\\u a — v°\\\ \ \v a \\i \\v a \\2 

It follows from (6.87), (6.88), and (6.79) that 

^_ap(w,w) + v\\Apw\\ 2 

< c(l + \\v°\\l + \\v a \\l)ap(w,w) (6.89) 
i 

+ccq \\u a — v°\\ 2 \\v a \\ 2 + ccoa||f a ||2 

Consequently, we can get 



f 1 

sup \\w(t)\\j + / \\w(t)\\ 2 dt 

0<t<Ti JO 

< f 'e^/o^d+INIIl+IKUDcit [ CCoa ||^(t)||2 +cc |||^_ v O||2|| v a|| 2 ]^ 

T 

< cia + ci f 1 \\u a -v \\l\\v a \ 
Jo 



< ci a + ci ^ 1 \\u a - v°\\idtj 2 1 \\v a \f 2 dt * 

Ti \ \ 



< Cl a + Cl (|[z;"||i + ||^ ||?) ( / |K -v°\\jdt ) < Cl a 



where we have used (6.75). Thus (6.76) holds, and the proposition is proved. 
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Remark 6.6 It is not clear to us whether the stronger estimate as (6.73) holds for u a 
under the assumptions in theorem 6. 3. . However, under the additional assumption that 
the strong solution v° £ L°°([0, T], H 2 ), there holds also 



r T i 

sup \\u a - v°\\l + / \\u a -v°\\ldt <ca2 (6.90) 

<t<T Jo 



ll« a -«°lli + / 

0<t<T 

This follows from 



a p {u a - v°, u a - v°) + a{PA(u a - v°),PA(u a - v )) 
= ap(v a - v°,u a - v°) + a(-PAv°, PA(u a - v )) 

(due to (6.83)) and (6.73). 

7 Concluding Remarks 

We conclude this paper with a few remarks on related issues. 

Remark 7.1 In exact same way, we can study the boundary value problem of LNS-a 
with NSB: 

d t v - Av + V x v x u + Vp = in 9, (7.1) 

V • v = in n (7.2) 
u — aAu + Vp = v in $7 (7.3) 

V • u = in ft (7.4) 
v ■ n = 0, 2(S(v)n) T = -jv T on dn (7.5) 
u ■ n = 0, 2{S{u)n) T = —ju T on <9fi (7.6) 

The functional setting is similar to that of (4- 1)-(4-6), and all the results stated in 
section 4-6 are also valid. 

Remark 7.2 The non-homogenous boundary value problems of LNS-a with VSB: 

d t v - Av + V x v x u + Vp = in Q. (7.7) 

V • v = in n (7.8) 
u — aAu + Vp = v in ft (7.9) 

V • u = in n (7.10) 
v ■ n = 0, n x V x v = (3v + b on dfl (7-11) 
u ■ n = 0, n x V x u = (3v + b on dtt (7.12) 

can also be considered by using a homogenous method to reduce it into 

d t v - Av + V x v x u + Vp = £ in n (7.13) 

V • v = in n (7.14) 
u — aAu + Vp = v + r] a in Q (7-15) 

V • u = in n (7.16) 
v . n = o, n x V x v = f3v on 90 (7.17) 
u • n = 0, n x V x u = f3u on dtt (7.18) 
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for some £ and n a as was done for the steady homogenous case in Section 3. Similarly, 
the non-homogenous boundary value problems for LNS-a with NSB may be established 
too. 



Remark 7.3 In the functional settings, the parameters associated with the velocity v 
and the filter u can be different, and different type boundary conditions, VSB or NSB, 
may be also allowed. However, in this case the analysis in the global existence and 
the vanishing a limit seems very difficult since 15.^7] ) does not hold, there are some 
boundary terms arising, and the energy estimate in \5.16\) depends on a. Yet the local 
well-posedness theory can be established by the method discussed in this paper. 

Remark 7.4 Our approaches works also for other a models. For instance, one can 
consider the following Leray a model: 

d t v- Av + u- Vv + Vp = in ft (7.19) 

V • v = in ft (7.20) 
u - aAu + Vp = v in ft (7.21) 

V • u = in ft (7.22) 
vn = 0, n x (V x v) = pv on 9ft (7.23) 
u ■ n = 0, 2(S(u)n) T = -ju T on 9ft (7.24) 

which allowed different boundary conditions between the velocity v and the filter u. In 
fact, this model is easier to analyze than the LNS-a since it has the following energy 
equation 

-f-IMI 2 + 2aR(v,v) = (7.25) 
dt 

which yields the global existence directly, and the corresponding convergence result is 
better both in v a and u a than that in theorem 5.1.. 
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